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ABSTRACT 

We explore the possibility of subsequent star formation after a first star forms in a Pop III object, 
by focusing on the radiation hydrodynamic (RHD) feedback brought by ionizing photons as well as 
II2 dissociating photons. For the purpose, we perform three-dimensional RHD simulations, where the 
radiative transfer of ionizing photons and II2 dissociating photons from a first star is self-consistently 
coupled with hydrodynamics based on a smoothed particle hydrodynamics method. As a result, it is 
shown that density peaks above a threshold density can keep collapsing owing to the shielding of II2 
dissociating radiation by an H2 shell formed ahead of a D-type ionization front. But, below the threshold 
density, an M-type ionization front accompanied by a shock propagates, and density peaks are radiation 
hydrodynamically evaporated by the shock. The threshold density is dependent on the distance from 
a source star, which is ~ lO^cm"^ for the source distance of 30pc. Taking into consideration that the 
extent of a Pop III object is ^ lOOpc and density peaks within it have the density of 10^~^cm~^, it is 
concluded that the secondary star formation is allowed in the broad regions in a Pop III object. 

Subject headings: theory: early universe — galaxies: formation — radiative transfer — molecular 
processes — hydrodynamics 



1. INTRODUCTION 

In the last decade, the formation of Population III (here- 
after Pop III) stars has been explored extensively (Bromm, 
Coppi & Larson 1999, 2002; Nakamura & Umemura 1999, 
2001; Abel, Bryan, Norman 2000, 2002; Yoshida 2006). 
The Pop III stars can significantly affect the reioniza- 
tion history in the universe (Cen 2003; Ciardi Ferrara, 
& White 2003; Wyithe & Loeb 2004; Somerville & Livio 
2003; Sokasian et al. 2004; Murakami et al. 2005), and the 
metal enrichment of intergalactic medium (Nakamura & 
Umemura 2001; Scannapieco, Ferrara, & Madau 2002; Ri- 
cotti & Ostriker 2004). In previous analyses (Tegmark et 
al. 1997; Fuller & Couchman 2000; Yoshida et al. 2003), 
it is shown that Pop III objects have the halo mass of 
^ lO^M© and the extent of ^ lOOpc. In Pop III objects, 
density peaks collapse owing to II2 cooling, forming cloud 
cores with the density of 10^~^cm~^ (Bromm, Coppi & 
Larson 2002). A highest peak in the halo collapses earlier 
to form a first Pop III star. Hence, subsequently collaps- 
ing cores are affected by the radiative feedback by the first 
star. However, the range of the feedback and the possibil- 
ity of secondary star formation in a Pop HI object are still 
under debate. 

If a first star is distant by more than Ipc, dense cores 
are readily self-shielded from the ultraviolet (UV) radia- 
tion (Tajiri & Umemura 1998; Kitayama et al. 2001; Susa 
& Umemura 2004a, b; Kitayama et al. 2004; Dijkstra et al. 
2004; Alvarez et al. 2006). Thus, the photoevaporation by 
UV heating is unlikely to work devastatingly. However, 
the photodissociating radiation (11. 18-13. 6eV) of hydro- 



gen molecules (H2) in Lyman- Werner band (LW band) can 
preclude the core from collapsing, since H2 are the dom- 
inant coolant to enable the core collapse (Haiman et al. 
1997; Omukai & Nishi 1999; Haiman et al. 2000; Glover 
& Brand 2001; Machacek et al. 2001). Hence, this may 
lead to momentous negative feedback. 

On the other hand, ionizing radiation for hydrogen 
(>13.6eV) drives an ionization front (I- front), which prop- 
agates in a collapsing core. The enhanced fraction of elec- 
trons promotes H2 formation (Shapiro & Kang 1987; Kang 
& Shapiro 1992; Susa et al. 1998; Oh & Haiman 2002). In 
particular, the mild ionization ahead of the I- front can gen- 
erate an H2 shell, which potentially shields H2 dissociating 
photons (Ricotti, Gnedin, & Shull 2001). This mechanism 
is likely to work positively to form Pop HI stars. In prac- 
tice, the propagation of I-front is complex. When UV irra- 
diates a dense core, the I-front changes from R-type on the 
surface to D-type inside the core. The transition occurs 
via an intermediate type (M-type), which is accompanied 
with the generation of shock (Kahn 1954). The shock can 
affect significantly the collapse of the core. This is a to- 
tally radiation hydrodynamic (RHD) process, which has 
not been hitherto explored in detail as feedback by first 
stars. 

In this Letter, we scrutinize the radiation hydrodynamic 
(RHD) feedback by a first star through the propagation of 
I-front in a dense core. For the purpose, we solve the 
three-dimensional RHD, where the radiative transfer of 
ionizing photons as well as H2 dissociating photons is self- 
consistently coupled with hydrodynamics. In §2, the key 



physical processes associated with I-front propagation are 
overviewed. The setup of simulation is described in §3, 
and numerical results are presented in §4. §5 is devoted to 
the conclusions. 



However, dissociating radiation could be shielded effec- 
tively if an H2 shell forms ahead of I-front. This process 
is tightly coupled with the propagation of I-front. 



2. BASIC PHYSICS 

If a collapsing core is irradiated by an ionizing source lo- 
cated at the distance D^ the propagation speed of I-front 
in the core is given by 

-IF = 21 km s ^.^^^^j^^ (^^3---3 

(1) 

where Nion is the number of ionizing photons per unit time 
and ncore is the number density in the cloud core. The 
sound speed in the regions cooling by II2 is ai ^ 1km s~^, 
while that in the ionized regions is a2 ^ 10km s~^. Hence, 
viF > 2a2 for low density cores or a nearby ionizing source, 
and therefore I-front becomes R-type. Supposing the cloud 
core size rcore is on the order of ai yTr/Gp^^, the propa- 
gation time of R-type front across the core satisfies 
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This means that an R-type front sweeps the core before 
the core collapses in a free-fall time. Thus, the core is 
likely to undergo photoevaporation. On the other hand, 
if the density of cloud core is high enough and the source 
distance is large, then vif < a1/2a2 and a D-type I-front 
emerges. The propagation time of D-type front across the 
core is 
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Thus, the core can collapse before the I-front sweeps the 
core. 

The above arguments are based on the assumption 
that the ionizing photon flux does not change during the 
propagation of I-front. However, the core could be self- 
shielded from the ionizing radiation if N{i:r1^^^/A'KD'^) < 
47rr^ore^core^B/3, whcrc q^b is the recombination coeffi- 
cient to all excited levels of hydrogen. The critical density 
for self-shielding is given by 
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If ncore > ^shield 7 the iouiziug photou flux diminishes sig- 
niflcantly during the I-front propagation. Hence, even if 
the I-front is R-type on the surface of cloud core, the front 
can change to M-type accompanied with shock, and even- 
tually to D-type inside the core (Kahn 1954). 

In comparison to ionizing photons, H2 dissociating (LW 
band) radiation is less shielded (Draine & Bertoldi 1996). 
The self-shielding of LW band flux (Flw) is expressed by 

Flw = ^LW,o/sh (^H2,14) (5) 

where Flw,o is the incident flux, A^H2,14 = A/'Hs/lO^^cm"^, 
and 

/sh(x) = I ^13/4 l%\ (6) 



3. SETUP OF SIMULATION 

We perform radiation hydrodynamic simulations using 
a novel radiation transfer solver based on smoothed par- 
ticle hydrodynamics (SPH) (Susa 2006). We suppose a 
primordial gas cloud collapsing in a run-away fashion. The 
chemical compositions are initially assumed to be the cos- 
mological residual value (Galli & Palla 1998). The mass 
of cloud is Mb = 8.3 x lO^M© in baryonic mass. When 
the central density of cloud exceeds a certain value, non, 
we ignite a I2OM0 Pop HI star (TVion = 1.3 x lO^^s"^), 
which is located dX D = 20pc for flducial models. We also 
explore the dependence on the distance of a flrst star by 
changing the relative location between a source star and 
a collapsing cloud. The luminosity and the effective tem- 
perature of the source star are taken from Baraffe, Heger 
& Woosely (2001). The ffducial models are 

Model A — non = 3 x lO^cm"^ with H2 dissociating 
photons but without ionizing photons 

Model B — non = 3 x lO^cm"^ with ionizing and H2 
dissociating photons 

Model C — same as Model B but for non = 3 x lO^cm"^ 
The physical simulation time after the ignition of the 
source star is 4Myr except 1.54Myr for model B, since 
the central part collapses below the resolution limit. In 
Model A, it is assumed that only LW band photons es- 
cape from the neighbor of a source star, whereas ionizing 
photons do not because of large opacity. This model also 
can be regarded as the reference to Model B. 



4. RESULTS 

4.1. Failed Collapse (Model A) 

In Fig. 1, the time evolution of density profiles along 
the axis of symmetry is shown for three models. In this 
figure, the red curve represents the distribution at IMyr, 
at which epoch the spatial distributions of physical quanti- 
ties are shown in Fig. 2. For Model A, as shown in the left 
panel of Fig. 1, the collapse of central regions stops vir- 
tually and forms a quasi- hydrostatic core between ~lMyr 
and ^4Myr. The failure of collapse is caused by the pho- 
todissociation of H2 by LW band photons. In upper panels 
of Fig. 2, the spatial distribution of H2 fraction on a slice 
along the symmetry axis is shown by colored dots at the 
positions of SPH particles. In lower panels of Fig. 2, phys- 
ical quantities along the symmetry axis are shown at this 
epoch. In Model A (left panel), the temperature is several 
lO^K in the envelope and ~ lO^K in the central regions. 
The H2 column density is lower than lO^^cm"^ in the en- 
velope and therefore H2 is highly photodissociated by LW 
band radiation. In the central regions, H2 column density 
exceeds lO^'^cm"^ and therefore H2 fraction is raised up 
by the self-shielding of LW band photons according to (5). 
However, this level of shielding is not enough to allow the 
cloud to keep collapsing by H2 cooling. Eventually, the 
collapse is halted by the thermal pressure. The numerical 
results of Model A show that gravitational instability is 
hindered by permeated H2 dissociating photons. 



4.2. H2 Shielded Collapse (Model B) 

In Model B, ionizing photons are included for the same 
model parameter as Model A. In this model, I-front prop- 
agates into the cloud. Although the I-front is R-type far 
from the cloud center, it changes to M-type and a shock is 
generated as shown by a peak at O.SMyr in Fig. 1 (middle 
panel). Around IMyr, the I-front changes to D-type, and 
therefore the core collapse proceeds faster than the prop- 
agation of the I-front, as argued in §2. The density distri- 
bution is highly changed by photoionization, as shown by 
the angelfish-shaped distribution in Fig. 2. In particular, 
it is worth noting that an II2 shell forms ahead of the I- 
front due to the enhanced ionization fraction. As a result, 
the H2 column density is steeply raised up. This H2 shell 
effectively shields LW band radiation from a source star, 
and the resultant H2 fraction {yn^) becomes larger by an 
order of magnitude than that of Model A. Eventually, due 
to H2 cooling, the cloud core can continue to collapse to 
the level of nn > lO^cm"^, as shown in Fig. 1 (middle). 
Since the only difference between Models A and B is the 
presence of ionizing radiation, we can conclude that pho- 
toionization can restrain the negative feedback effect by 
II2 photodissociation. 

4.3. Shock- driven Evaporation (Model C) 

Model C is the lower density version of Model B, where 
non is ten times smaller than that of Model B. In this case, 
quite similarly to Model B, the I-front becomes to M-type 
at O.SMyr (Fig. 1, right panel). As shown in Fig. 2, there 
forms an II2 shell that shields LW band radiation from a 
source star. But, before the core collapses in the free-fall 
time, the shock accompanied with M-type I-front sweeps 
up the central regions at ^ IMyr. Eventually, the shock 
blows out the collapsing core. Hence, it is concluded that 
if an M-type I-front passes through the cloud core, the 
whole cloud is evaporated radiation hydrodynamically by 
the shock. 

4.4. Dependence on Position of a Source Star 

As shown above, the criterion of radiation hydrody- 
namic evaporation of a collapsing core is whether an M- 
type I-front sweeps up the core. By changing the position 
of a source star, we can obtain the threshold density for 



the evaporation depending on the source distance. As a 
result of numerical simulations, it is shown that the thresh- 
old density is ~ lO^cm"^ for D = 20pc, ~ lO^cm"^ for 
D = 30pc, and ^ 10cm~^ for D = 50pc. Thus, cloud cores 
with ncore ^ lO^cm"^ can collapse at I) > 30pc. Further 
details of the analysis on the parameter dependence will 
be described in the forth-coming full paper. 

5. CONCLUSIONS 

We have performed three-dimensional radiation hydro- 
dynamic simulations to scrutinize the feedback by a first 
star in a Pop III object. As a result, it has been found 
that a collapsing core is evaporated by a shock if an M- 
type I-front sweeps the core. In order for a collapsing 
core to evade the radiation hydrodynamic evaporation, 
the core density should exceed a threshold density that 
depends on the distance from a source star. Above the 
threshold density, the I-front is quickly changed to the D- 
type, and an H2 shell forms ahead of the I-front, which 
effectively shields H2 dissociating radiation from a source 
star. Eventually, the core can keep collapsing owing to H2 
cooling. The present numerical study has shown that if a 
source star is distant by more than 30pc, then a collapsing 
core denser than ^ lO^cm"^ is absolved from evaporation. 
Taking into account that the density of cloud cores in a 
Pop III object is expected to be 10^~^cm~^, stars can form 
in the regions of > 30 pc even after a first star forms. Since 
the extent of a Pop III object is ^ lOOpc, it is concluded 
that the subsequent star formation is allowed in the broad 
regions in a Pop III object. 
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Fig. 1. — The time evolution of density distributions along the symmetry axis for three models. Four curves correspond to t = 0, 0.5Myr, 
IMyr, and 4Myr, for models A and C, whereas the final plot for model B corresponds to t = 1.5Myr. The red curves denote the profiles at 
IMyr, at which epoch the detailed structure is shown in Fig. 2. 
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Fig. 2. — Spatial distributions of physical quantities for three models at 1 Myr. Colored dots on upper three panels represent the distributions 
of SPH particles on a slice which includes the symmetry axis (dashed line). The colors of particles denote the H2 fraction, whose legend is 
shown on the upper left. The yellow star represents the position of a source Pop III star. Lower three panels show various physical quantities 
along the symmetry axis as TlK], nH[cm~^], ym^ 2/H25 ^"^^ A^H2, 14(^2 column density from the source star in units of lO-'^^cm"^). Horizontal 
axis shows the distance from the source star. Arrows show the correspondence of coordinates between the SPH distributions and the horizontal 
axis. 



